Temporal distortion of annual modulation at low recoil energies 
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We show that the main features of the annual modulation of the signal expected in a WIMP direct 
detection experiment, i.e. its sinusoidal dependence with time, the occurrence of its maxima and 
minima during the year and (under some circumstances) even the one-year period, may be affected 
by relaxing the isothermal sphere hypothesis in the description of the WIMP velocity phase space. 
The most relevant effect is a distortion of the time-behaviour at low recoil energies for anisotropic 
galactic halos. While some of these effects turn out to be relevant at recoil energies below the 
current detector thresholds, some others could already be measurable, although some degree of 
tuning between the WIMP mass and the experimental parameters would be required. Either the 
observation or non-observation of these effects could provide clues on the phase space distribution 
of our galactic halo. 
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INTRODUCTION 

It is widely believed, as suggested by a host of indepen- 
dent cosmological and astrophysical observations, that 
the most part of the matter in the Universe is not visible, 
revealing its existence only through gravitational effects. 
In particular, data on the rotational curves of galaxies 
indicate that the galactic visible parts are surrounded by 
approximately-spherical dark halos which extend up to 
several times the size of the luminous components. 

The best candidates to provide dark matter in galaxies 
are Weakly Interacting Massive Particles (WIMP). Sev- 
eral WIMP direct-detection experiments are operating 
, with the goal of measuring the nuclear recoil energy 
(in the KeV range) expected to be deposited in solid, 
liquid or gaseous targets by the scattering of the non- 
relativistic dark halo WIMPs. Unfortunately, expected 
rates are small and the exponential decay of the WIMP 
recoil-spectrum resembles that of the background at low 
energies. However, a specific signature can be exploited 
in order to disentangle a WIMP signal from the back- 
ground: the annual modulation of the rate Q . This ef- 
fect, expected to be of the order of a few per cent, is 
induced by the rotation of the Earth around the Sun. 
Due to its smallness, the annual modulation signature 
requires large-mass detectors with high statistics in or- 
der to overcome background fluctuations and be unam- 
biguously detected. The annual modulation effect has 
been experimentally investigated by the DAMA Collab- 
oration, which has indeed reported a positive evidence 
by using a 100 kg sodium iodide detector 

One of the most important sources of uncertainty in 
the calculation of WIMP direct detection rates is the 
modeling of the velocity distribution function (DF) of 
the particles populating the dark halo. In the literature, 
a simple isothermal sphere model is usually adopted, i.e. 



a WIMP gas described by an isotropic Maxwcllian with 
r.m.s velocity of the order of 300 Km s _1 . This leads to 
a sinusoidal time-dependence of the expected signal with 
maximum (or minimum) around June 2 nd , i.e. with the 
same (or opposite) phase as the relative velocity between 
the Earth and the halo rest frame. 

However, the actual form of the WIMP velocity DF is 
unknown, and many different models, alternative to the 
isothermal sphere, are compatible with observations [jj. 
The goal of the present Letter is to show that the main 
features of the annual modulation effect (the sinusoidal 
dependence with time, the occurrence of maxima and 
minima during the year and, under some circumstances, 
the one-year period) may be affected by anisotropics in 
the velocity DF. The most relevant effect is a distortion 
of the sinusoidal time-behaviour at low recoil energies. 
These energies, though below the current detector thresh- 
olds, might be reached in the future. The observation of 
the effects discussed in this Letter could provide infor- 
mations on the phase space distribution of our galactic 
halo, especially on the degree of its anisotropy. 



THE ANNUAL MODULATION EFFECT 

Due to the rotation of the disk around the galactic 
center, the solar system moves through the WIMP halo, 
assumed to be at rest in the galactic rest frame. In the 
following, we will assume a right-handed system of or- 
thogonal coordinates: the x axis in the galactic plane, 
pointing radially outward; the y axis in the galactic plane, 
pointing in the direction of the disk rotation; the z axis 
directed upward, perpendicular to the galactic plane. No- 
tice that our system differs from standard "galactic co- 
ordinates" by the different choice of the x axis, which for 
us is directed outward. 
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The relative velocity between the WIMP halo and the 
detector is given by the Earth velocity v E , as seen in the 
galactic rest frame. It is the sum of three components: 
the galactic rotational velocity = (0,v ,0) Km s _1 
(we will assume: vq = 220 Km s _1 ), the Sun proper 
motion v s = (-9, 12, 7) Km s _1 [3| and the Earth orbital 
motion u s (/j)[4|]: 



+ v% + u E (\)cosf3 x cos[Lo(t-t x )], (1) 
+ v% + u E (\)cos(3 y cos[u(t-t y )}, (2) 
v s z + w £ (A)cos/3 z cos[w(t-^)], (3) 
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where A is the ecliptic longitude, which is function of 
time. We can express A as Q: A = L + 1°. 915 sing + 
0°.020sin2 5 where L = 280°.460 + 0°. 9856474 t and 
g = 357°. 528 + 0°.9856003 t and t denotes the time ex- 
pressed in days relative to UT noon on December 31. In 
Eqs.JTHSJ u E (X) = (u B )[l-esin(A- A )] is the modulus 
of the Earth rotational velocity, which slightly changes 
with time due to the small ellipticity e of the Earth orbit 
u E ) = 29.79 Km s" 1 , e = 0.016722 and A = 13° ± 1°) 
In Eqs.mOJ) the Pi denote the ecliptic latitudes 
and the ti are the phases of the three velocity compo- 
nents: f3 x = 174°. 4697, P y = 59°. 575, p z = 29°. 812 and 
t x = 76.1 day, t y = 156.3 day, t z — 352.4 day. The an- 
gular velocity has a period of 1 year, and is given by: 
uj = 27r/(365 days). With the numbers given above, 
the modulus of the Earth velocity changes in time as: 
v E = \v E \ ^ 233. 5 + 14. 4cos[u(t-t )] (inKms" 1 ), where 
t — 152 days, i.e. June 2 nd . Notice the slight offset be- 
tween to aud t y , due to the composition of the velocity 
components. As a good approximation, ve is usually 
taken as: ve = (vo + Vy ) + (u E ) cos P y cos[w(t — i )]. 

The direct detection differential rate dR/dEn is pro- 
portional to the integral: 



dw 



/es(w) 



(4) 



w>v n 



where /es and w are the WIMP velocity DF and the 
WIMP velocity, respectively, in the Earth's rest frame; 
v m in is the minimum value of w for a given recoil energy 
En, WIMP mass mw and nuclear target mass mjv and 
is given by: v min = y/ E R J (2m N ){m w + m N )/m w . 

By indicating with f(v) and v the WIMP velocity DF 
and the WIMP velocity in the galactic reference frame, 
the following transformations hold: 

v — ► w = v — v E (t), (5) 
m - f ES (w) = f(n! + v E (t)), (6) 




FIG. 1: Time dependence of X(v m i n ) for an isotropic isother- 
mal sphere. The different curves refer to values of v m i n rang- 
ing from (lower) to 400 (upper) Km s _1 , in steps of 20 Km 
s _1 (the dashed curves correspond to w m i n = 100, 200 and 300 
Km s _1 ). The vertical dashed line denotes t = 152 days. 



THE ISOTROPIC ISOTHERMAL SPHERE 

In the case of the isothermal sphere model the DF is 
given by a truncated isotropic Maxwcllian which depends 
only on v = \v\: 



f(v) = ./Vexp 



v 

2^2 



(7) 



where a is the WIMP r.m.s. velocity, given by: a 2 = 
3«o/2. It is clear that, through the change of refer- 
ence frame of Eq. 2"(i> m in) depends on time only 
through ve- Since the relative change of ve during the 
year is of the order of a few per cent, we can approxi- 
mate X(v m i n ) with its first-order expansion in the small 
parameter e = 8ve/ve, around its mean value To(v m i n ): 
Z(vmin) — 2"o( u min) oc cos[u>(t — to)]. The well known re- 
sult is then obtained that the WIMP rate has a sinusoidal 
time dependence with the same phase (to — June 2 nd ) as 
ve, for all values of u m in- This is shown in Fig^ where 
X(u m i n ) is plotted as a function of time for various values 
of u mi „. 



ANISOTROPIC MODELS 

The simplest generalization of the isothermal sphere 
model is given by a triaxial system described by a mul- 
tivariate gaussian: 



which imply that T(v m i n ), and so dR/dEn, develops a 
time dependence induced by VE(t). 



f(v) = N exp 



2ai 



2a? 



(8) 
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FIG. 2: The same as in Fig. for an anisotropic model with 
A12 = tJx/o-y — 0.2 and A32 = o- z /a y — 0.8. 

where TV is the normalization constant. For v csc — > 00, 
then N = [(2Tral)(2ira 2 )(2Tra 2 )}- 1 / 2 . The usual isother- 
mal sphere is the spherical limit of Eq. JSJ) , obtained with: 
a x = Cy = a z = In order to discuss the effect 
of anisotropy at fixed WIMP mean kinetic energy, we 
will fix a 2 = a 2 . + a 2 + a 2 as in the isothermal case 
(a 2 = 3v 2 /2) and discuss our results in terms of the two 
independent parameters: A12 = cr x /(Ty and A32 = cr z /a y . 

At variance with the isothermal sphere, now I(v m in) 
depends in general on all the three components of ve, 
and not simply on ve- We can write: 

f ES (w) = Ne X p[-(r + rj) 2 ] , (9) 

where we have defined the reduced (dimensionless) vari- 
ables: n = w t /(V2<r t ) and rji = vf/(y/2ai) = rft + 
Srji cos[u(t — ti)] (i = x,y,z). In the isotropic case 
(isothermal sphere) one has rj x z <C 1, 77^ ~ 1 and 
Si]i -c 1. The presence of the order-one parameter 77° 
and of the small oscillation amplitudes 5rji allows the 
Taylor-expansion of I(v min ) in terms of ej = drji/rpy pa- 
rameters. A straightforward calculation shows that the 
conclusions of the previous Section are recovered. 

On the contrary, allowing anisotropics such that A 12 < 
1 and/or A32 < 1, the values of the parameters r\ XiZ and 
Sr] x .z are enhanced, and the time-dependence of and 
vf in X(u m in) may become important. An example of this 
situation is shown in Fig. [3 where, the time evolution 
of l(v m i n ) is plotted for A12 = 0.2 and A12 = 0.8 (i.e.: 
(a x ,a y ,a z ) = (42, 208, 166) Km s _1 .) This choice of A i2 , 
A3 2 , which refers to a tangential anisotropy, corresponds 
to triaxial models discussed, for instance, in Refs. 
A distortion of the curves of Fig. as compared to 
the familiar sinusoidal time-dependence, appears: this 
effect may be explained by the fact that now 5rj x ~ 1, 
a Taylor-expansion of the type used in the isothermal 



FIG. 3: The same as in Fig. for an anisotropic model with 
A12 = <Jx/o-y — 10 and A32 = o z ja v = 3. The curves refer to: 
u min = (lower), 180, 190, 200 (dashed), 210, 220, 230, 240, 
250 (upper) Km s" 1 . 



sphere case breaks down and a full numerical calculation 
of the integral of Eq.QJ is required. The final result is not 
sinusoidal. This peculiar behaviour is more pronounced 
at low values of v m i n (i.e. low recoil energies) namely for 
^min < 80 Km s" 1 , for which the distortion is strong and 
the maxima (in absolute value) of the rate are shifted as 
compared to the standard case 0- For larger values of 
Vmin the distortion is less pronounced, and it dies away 
when v mnl > 230 Km s _1 . This may be explained by the 
fact that, as v m i n grows, the integral of Eq.Q becomes 
less sensitive to the parameter rj x since it gets increasingly 
dominated by WIMPs with velocities along the y axis, 
which is the one along which the boost due to the galactic 
rotation is directed. We notice that for values of V m i B 
around (100—200) Km s _1 a distortion is present, but the 
amplitude of the modulation is suppressed and therefore 
difficult to detect. 

As a second example, in FigO we plot 2(v m - nl ) as a 
function of time for the case A12 = 10, A 32 = 3 (i.e.: 
(a x ,(Ty,a z ) — (257, 26, 77) Km s _1 . This situation is rep- 
resentative of a radial anisotropy. In this case, we have 
further enhanced the contribution of Srjy over Sr] x ^ Z7 so 
that one should expect to draw the same conclusions as 
in the isothermal sphere case, with the usual sinusoidal 
time-dependence of the rate and a phase close to t y ~ t . 
This is indeed the case, except for a narrow interval of 
values of v m - nl around v mm ~ 210 Km s^ 1 . In this range, 
2(vmin) develops two maxima because, for that particular 
choice of u m i n , there is an exact cancellation in the first 
term of the expansion in e y = Sr/y/rjy, so that the term 
0(e 2 ) proportional to cos 2 [uj (t — t y )] sets in. This partic- 
ular cancellation of the first term in the Taylor expansion 
of T(v m i n ) happens also for the isothermal-sphere model, 
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FIG. 4: Values of v m i n as a function of the quenched nuclear 
recoil energies E ee = qEn, for: Nal scintillators 0, Ge ion- 
ization detectors Q] and Te bolometers For each panel, 
the different curves refer to WIMP masses of: 20 GeV (solid), 
50 GeV (dotted), 100 GeV (dashed), 200 GeV (long-dashed) 
and 1 TeV (dot-dashed). The dashed vertical lines denote 
the current energy thresholds. 



but in that case the size of the quadratic term is strongly 
suppressed because e y is much smaller. In a WIMP di- 
rect detection experiment this effect would show up in a 
very peculiar way: a halving of the modulation period 
of the rate in a narrow range of recoil energies. Note, 
however, that in order to have some realistic chance to 
detect this effect, it should show up in one of the ex- 
perimental energy bins just above threshold, where the 
highest signal/background ratio is usually attained. 

In order to establish a link between our discussion and 
WIMP direct detection experiments 0, in Fig. 0|we plot 
v m in as a function of the quenched nuclear recoil energy 
E ee = qEf> (q is the quenching factor) for the target 
nuclei: Na, I, Ge, Te and for different WIMP masses. 
The vertical dashed lines show current energy thresholds 
achieved by each type of detector. Fig. 0]shows that val- 
ues of f m ; n < 80 Km s — 1 , i.e. sufficiently low to observe 
a sizeable distortion effect as the one discussed for tan- 
gential anisotropy, correspond to WIMP recoil energies 
below the threshold of present direct detection experi- 
ments, and that the effect would be more easily detected 
at higher WIMP masses. However, a foreseeable lower- 
ing of the threshold, down to 0.5-1 KeV, could be enough 
to observe the distortion. On the other hand, for radial 
anisotropy, we can conclude that the recoil energy cor- 
responding to a halving of the modulation period can 
actually coincide to the experimental thresholds within 
the reach of present-day detectors for 20 GeV < my/ < 
100 GeV, depending on the particular target nucleus. In 



this respect, wc note that the properties of the annual 
modulation effect observed by the DAMA/Nal experi- 
ment [j| (a one-year-period sinusoidal behaviour in the 
2-6 KeV energy bins implies that the DAMA/Nal 

experiment is already able to set constraint on strong 
radial anisotropies. 



CONCLUSIONS 

In the present Letter we have shown that the main 
features of the annual modulation of the signal of WIMP 
direct searches, i.e. the sinusoidal dependence of the rate 
with time, the position of its maxima and minima during 
the year and even the period, may be affected by relax- 
ing the isothermal sphere hypothesis in the description of 
the WIMP velocity phase space. We have considered a 
multivariate gaussian and found that different situations 
may occur, depending on the pattern of anisotropy: tan- 
gential anisotropies induce a departure at low energies 
from the usual sinusoidal time-dependence, along with a 
shift in the position of the maximum of the signal during 
the year, while radial anisotropics may produce a halving 
of the modulation period in a particular energy bin. The 
former effect turns out to be relevant at low recoil ener- 
gies, actually below the threshold of present-day experi- 
ments, while the latter should be already within the reach 
of current detectors. In particular, the properties of the 
annual modulation effect observed by the DAMA/Nal 
experiment Q may already indicate that strong radial 
anisotropies are excluded. 

The effects discussed in this Letter may be used in the 
future to provide a direct way of measuring (or setting 
limits on) the degree of anisotropy of the galactic halo. 
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